Rotation numbers for Jacobi matrices with matrix entries 
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Abstract 

A selfadjoined block tridiagonal matrix with positive definite blocks on the off-diagonals 
is by definition a Jacobi matrix with matrix entries. Transfer matrix techniques are extended 
in order to develop a rotation number calculation for its eigenvalues. This is a matricial gen- 
eralization of the oscillation theorem for the discrete analogues of Sturm-Liouville operators. 
The three universality classes of time reversal invariance are dealt with by implementing 
the corresponding symmetries. For Jacobi matrices with random matrix entries, this leads 
to a formula for the integrated density of states which can be calculated perturbatively in 
the coupling constant of the randomness with an optimal control on the error terms. 
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1 Introduction 

This article is about matrices of the type 

/ V 1 T 2 

T 2 V 2 

T 3 



H 



N 



V 
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V 3 



\ 



Vjv-i Tjsi 
T N V N J 
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where V n = V* are selfadjoint complex Lx L matrices and T n are positive definite complex Lx L 
matrices (Note added Feb. 2008: invertibility of T n is sufficient with minor modifications). With 
the convention T\ = 1 and for a complex energy E G C, introduce the transfer matrices 
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Theorem 1 Let E G K and N > 2. 

(i) U§ is well-defined, namely the appearing inverse exists. 

(ii) U§ is a unitary matrix which is real analytic in E. 

(iii) The real eigenphases 6^ v I = 1, . . . , L, of U§ can be chosen (at level crossings) to be 
analytic in E and such that — > as E — > — oo and — > 2ttN as E — > oo. 

(iv) E is an eigenvalue of H of multiplicity m if and only if 9f Jl = ir mod 2tt for exactly m 
of the indices I = 1, . . . , L. 

(v) The matrix S§ = - (U§)*8eU§ is positive definite. Each 9f }l is an increasing function of E. 

(vi) If H N is real, the unitary U§ is symmetric and the positive matrix S§- is real. 

(vii) Let L be even and let 

1 = ( i ~o ) G Mat ^ L x L,C) ^ 

with 4 blocks of size ~ X ~ . Suppose that H N is self-dual, namely the entries are self-dual: 

PT l n I = T n , rV*I = V n , n = l,...,N. 

Then and S§ are also self-dual (equivalently, IU§ and IS§ are skew-symmetric) . 

Items (i), (ii), (iv), (vi) and (vii) result directly from the mathematical set-up, while the 
analyticity statements of items (ii) and (iii) are based on elementary analytic perturbation theory 
[KatJ. The second part of (iii) follows from a homotopy argument and item (v), even though 
a consequence of a straight-forward calculation, is the main mathematical insight. It justifies 
the term rotation numbers for the eigenphases. In the strictly one-dimensional situation and for 
Sturm-Liouville operators instead of Jacobi matrices, the theorem has been known for almost 
two cenuries as the rotation number calculation or the Sturm-Liouville oscillation theorem [Weil 
IJMj . For matricial Sturm-Liouville operators, Bott [Bot] has proven results related to the above 
theorem. (The author learned of Bott's work once this article was finished, and believes that the 
techniques presented below allow to considerably simplify Bott's proof. A detailed treatment is 
under preparation.) For related work on linear Hamiltonian system let us refer to the review 
[FJNj . The discrete one-dimensional case and hence precisely the case L — 1 of Theorem [T] 
is also well-known (see e.g. |JSSj for a short proof). A rougher result was proven by Arnold 
[Arn2j . In the one- dimensional situation the variable 8§ l is called the Priifer phase. Therefore 
one may refer to the eigenphases 9fj l (or the unitaries U§ themselves) also as mult i- dimensional 
Priifer phases. The two supplementary symmetries considered in items (vi) and (vii) correspond 
to quantum-mechanical Hamiltonians H with time-reversal invariance describing systems with 
odd or even spin respectively |Meh] . This notion is empty in the one-dimensional situation where 
time-reversal invariance follows automatically from self-adjointness. 

Crucial ingredient of the proof is that ([3]) for real energies actually stems from the Mobius 
action of the symplectic transfer matrices (j5J) on the unitary matrices (Theorem [5]), which in turn 
are diffeomorphic to the Lagrangian Grassmannian via the stereographic projection (Theorem H]). 
As a function of real energy, U§ hence corresponds to a path of Lagrangian planes. If one 
defines a singular cycle in the unitary group as the set of unitaries with eigenvalue —1, then the 
intersections of the above path with this cycle turn out to be precisely at the eigenvalues of H N . 
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One new perspective opened by the Theorem [T] concerns Jacobi matrices with random matrix 
entries, describing e.g. finite volume approximations of the higher- dimensional Anderson model. 
In fact, the unitary, symmetric unitary and anti-symmetric unitary matrices form precisely the 
state spaces of Dyson's circular ensembles. They are furnished with unique invariant measures 
(Haar measures on the corresponding symmetric spaces). A good working hypothesis is hence 
that the random dynamical system induced by the action of the symplectic transfer matrices 
on the unitary matrices has an invariant measure (in the sense of Furstenberg [BL] ) which is 
invariant under the adequate action of the unitary group. This can only be true to lowest order 
in perturbation theory, under a hypothesis on the coupling of the randomness (which has to be 
checked for concrete models), and on a set of lower dimensional unitary matrices corresponding 
to the elliptic channels (in the sense of [SBj and Section HI). This then would justify the random 
phase approximation or maximal entropy Ansatz (here as equidistribution of Lagrangian planes) 
widely used in the physics community in the study of quasi-one-dimensional systems and in order 
to establish a link between random models like the Anderson Hamiltonian and invariant random 
matrix ensembles (e.g. [Bee]). Furthermore, let us consider H describing a physical system on 
a (i-dimensional cube, namely with L = iV^" 1 , and suppose d sufficiently large. Then a further 
working hypothesis is that the positive matrices jjS§ are distributed according to the Wishard 
ensemble of adequate symmetry (again on the elliptic subspaces and in the weak coupling limit). 
If both working hypothesis turn out to hold and U§ and Sfj- are asymptotically independent, 
namely randomly rotated w.r.t. each other (for large cubes), then Theorem [1] combined with 
a convolution argument shows that the level statistics of H N is asymptotically given by the 
Wigner-Dyson statistics for quantum systems without or with time-reversal invariance for odd 
or even spin, according to the symmetry of H . Roughly, Theorem [1] hence gives one possible 
way to materialize the heuristics given in the introduction to Chapter 9 of Mehta's book |Meh] . 
Numerical results supporting the above have been obtained in collaboration with R. Romer. 

As a first application of Theorem [1] and the techniques elaborated in its proof we develop 
in Section H] the lowest order perturbation theory for the integrated density of states (IDS) 
of a semi-infinite real Jacobi matrix with random matrix entries. More precisely, we suppose 
that the entries in ([T} are of the form V n = V(l + Xv n + 0(X 2 )) and T n = T(l + \ t n + 
0(X 2 )) where V,T, v n , t n are real symmetric matrices and T is positive definite. The v n ,t n are 
drawn independently and identically from a bounded ensemble (v a , t (J ) - e s according to a given 
distribution E CT , and furthermore the dependence on the coupling constant A > is real analytic 
and the error terms satisfy norm estimates. Associated to a random sequence u = (v n , t n ) n >i are 
random real Hamiltonians H N (u, A). The Anderson model on a strip is an example within this 
class of models. The number of eigenvalues of H N (u, A) smaller than a given energy E and per 
volume element NL is a self-averaging quantity in the limit N — > oo which converges to the IDS 
J\f\(E) (see Section H] for the formal definition). Let furthermore M\, a (E) denote the IDS of the 
translation invariant Hamiltonian with uo = (v a ,t a ) n >i. Finally let T E be the transfer matrix 
defined as in (j2J) from the unperturbed entries V and T. 

Theorem 2 Suppose that E e R is such that T E is diagonalizable and does not have anomalies, 
namely the rotation phases of the elliptic channels are incommensurate (cf. Section [4721 for the 
precise hypothesis). Then 

M X (E) = E a (Afx, a (E)) + 0(\ 2 ) . 
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The same result holds for random perturbations of arbitrary periodic operators. The fact 
that T E is not allowed to have any Jordan blocks means that E is not an internal band edge. 
Together with the anomalies they form a discrete set of excluded energies. The main point of 
Theorem [2] is not the calculation of the leading order term in A (which is indeed given by the 
most naive guess), but rather the control of the error term which is uniform in L as long as 
one stays uniformly bounded away from anomalies and internal band edges. The estimates in 
Section 14.21 also show how the error bound diverges as one approaches these energies. However, 
this part of the analysis is not optimized and there is a definite need for refinement in order 
to be able to study the thermodynamic limit of solid state physics models. The error bound 
is nevertheless optimal in the sense that the 0(X 2 ) contribution on the l.h.s. does depend on 
further details of the model. Similar as in the one-dimensional situation (L = 1), the IDS and 
the sum of the positive Lyapunov exponents are imaginary and real boundary values of a single 
Herglotz function |KSj . Hence we also develop a perturbation theory for the sum of the Lyapunov 
exponents, with a considerably better control on the error terms than in [SBJ where a particular 
case has been treated. 

This work is organized as follows. The next section recollects the tools from symplectic 
geometry used in the proof of Theorem [U In particular, it is shown that the Mobius action of 
the symplectic group on the unitary matrices is well-defined and furthermore some formulas for 
the calculation of the intersection number (Maslov index) are given. Section [3] provides the proof 
of Theorem [1] and then gives some supplementary results on Jacobi matrices with matrix entries 
and their spectra. Section H] contains the definition of the IDS for Jacobi matrices with random 
matrix entries and a formula for the associated averaged Lyapunov exponent. Then the proof of 
Theorem [2] is given. 

Acknowledgment: The author thanks Hajo Leschke, Demetris Pliakis and Robert Schrader for 
discussions on the matters of the paper. This work was supported by the DFG. 



2 Symplectic artillery 

Apart from the propositions in Section 12.81 which may be strictly speaking new, this section is 
probably known to the experts in symplectic geometry. But there does not seem to be reference 
with a treatment as compact and unified as the present one. The author's references were 
[Hual ISiel lArnlt ICH IKoSl IArn3j and he hereby excuses for not citing all the interesting works 
that he does not know of. The reader is warned that the complex Lagrangian planes and the 
complex symplectic group are defined with the adjoint rather than the transpose. This differs 
from standard references, but hopefully the reader will agree that it is natural in the present 
context. 

Let us introduce some notations. The following 2L x 2L matrices (matrices of this size are 
denoted by mathcal symbols in this work) are composed by 4 blocks of size L x L: 

'-(:-)■ '-(i-°0- '-Mi".! 1 )- *-(; - 

where in the last equation I is given by (jlj) and hence L is supposed to be even. J is called 
the symplectic form, C the Cayley transform and I the self-duality transform. The following 



4 



identities will be useful: 

CJC* = - Q, CJC* = - J, CXC* = -X. (5) 

i % % 

In order to deal with the symmetry of Theorem (TJ^vii), hence L even, some further notations 
are convenient. A matrix A G Mat(L x L, C) is call self-dual if I*A t I = A, and it is called 
self-conjugate if I*AI = A. As already indicated in Theorem DJvii), self-duality is closely linked 
to skew-symmetry, namely A is self-dual if and only if {I A) 1 = —I A. The sets of skew-symmetric 
and self-dual matrices are denoted by Skew(L, C) and Self(L, C) respectively. Moreover, for 
selfadjoint matrices A* = A the notions of self-duality and self-conjugacy coincide. 

2.1 Lagrangian Grassmannian 

The set of L-dimensional subspaces of the complex vector space C 2i is denoted by G^. The 
vectors of a basis of such a plane form the column vectors of a 2L x L matrix <3> which has rank 
L. Of course, a plane does not depend on the choice of the basis (and hence the explicit form of 
$). Consider the relation: $ ~ ^ <^ there exists c G G1(L, C) with $ = \l/c. The Grassmannian 
is then the set of equivalence classes w.r.t. ~: 

= {[$]„ | $ G Mat(2L x L,C) , rank($) — L } . 
A plane is called (complex hermitian) Lagrangian if = 0. Here A* = A 1 denotes 



transpose of the complex conjugate of a matrix A. If $ = 




where a and b are complex 



Lx L matrices, the latter condition means that a*b = b*a is selfadjoint. The (complex hermitian) 
Lagrangian Grassmannian is the set of Lagrangian planes: 

= {[$]„ | $ G Mat(2L x L, C) , rank($) = L , <$>*J<$> = } . (6) 

This is a real analytic manifold. It contains the submanifold Lf of real Lagrangian planes: 

= {[$]„ | $ G Mat(2L x L, C) , rank($) = L , <&*J<& = , = } . (7) 

Hence is a subset of characterized by a supplementary symmetry. That this coincides with 
the usual definition of the real Lagrangian Grassmannian is stated in Theorem [3] below. If L is 
even, then contains another submanifold characterized by another symmetry: 

L= = | $ G Mat(2L x L, C) , rank($) = L , <$>*J<$> = , $*J$ = } . (8) 

The notation appealing to the quaterions is justified by the following theorem, in which H L 
is considered as a vector space over C. Let A* B denote the transpose and quaternion conjugate 
(inversion of sign of all three imaginary parts) of a matrix A with quaternion entries. Similarly, 
A* R = A 1 for a matrix with real entries. 

Theorem 3 One has, with equality in the sense of diffeomorphic real analytic manifolds, 
= {[$]„ | $ G Mat(2L x L,R) , rank($) = L , = } , 

and 

= {[$]„ | $ G Mat(L x L,H) , rank($) = L , = } . 

The proof is postponed to the next section. 
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2.2 Stereographic projection 



Bott [Bot] showed that the complex Lagrangian Grassmannian L 1 ^ is homeomorphic to the unitary 
group U(L), a fact that was rediscovered in |KoSt IArn3j . Arnold |Arnlj used the fact that Lf 
can be identified with U(L)/ O(L). Indeed, a real Lagrangian plane can always be spanned by an 
orthonormal system, that is, be represented by $ satisfying $*$ = 1. This induces ($, J7"$) G 
SP(2L, R)n 0(2L) = U(L). Of course, various orthonormal systems obtained by orthogonal basis 
changes within the plane span the same Lagrangian plane. Hence Lf = U(L)/ O(L). Moreover, 
the symmetric space U(L) / 0{L) can be identified with the unitary symmetric matrices by sending 
a right equivalence class AO(L) G U(L)/ O(L) to AA l . As these facts will be crucial later on, 
let us give a detailed proof and some explicit formulas. 

The stereographic projection tc is defined on the subset 



G G L | (01) $ G GL(L,C)} 



by 



tt([$]. 



(10) $ ( (01) $) 1 



a b 



One readily checks that 7r([$]^) is independent of the representative. If [$]^ G L^ n G£ v , then 
7r([$]^) is selfadjoint. The fact that tt is not defined on all of is an unpleasant feature that 
can be circumvented by use of II defined by 



n([*] r 



Ti 



([CS], 



if [C$]„ G cr 



Theorem 4 (i) T/ie maj> LT : L 1 ^ — ►U(L) is a real analytic diffeomorphism. 

(ii) TTie map U : Lf — ►U(L) fl Sym(L, C) is a rea/ analytic diffeomorphism. 

(iii) Lei L 6e even. TTie maj> LT : Lf — » U(L) fl Self(L, C) is a rea/ analytic diffeomorphism. 



Proof. Let $ 



where a and b are L x L matrices satisfying a*b = b*a. One has 



rank(<3>) = rank($*$) = rank(a*a + b*b) 



rank((a + i b)*(a + 1 6)) = rank(a + 2&) = rank (a — ib) . 



(9) 



It follows that G G£ v so that it is in the domain of the stereographic projection n and 

hence LT is well-defined. Next let us show that the image is unitary. It follows from Q and a 
short calculation (or alternatively the first identity in (jSJ)) that 



Ch c L 



a,be GL(L,C) , a*a = b*b 



Hence if [C$] r 



G C L£, one has 



II([$]„)*II([$]„) = (b*)- 1 a*ab~ 1 



1 . 
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Moreover, II is continuous. One can directly check that the inverse of II is given by 



u-\u) 



K&r + i) 



As this is moreover real analytic, this proves (i). 

For the case (ii) of the real Lagrangian Grassmannian, the second identity of 
the supplementary symmetry in (jTJ) leads to 



(10) 



implies that 



C L, 



a,6 G GL(L, C) 



a a 



b*b , (ah 



~i\t 



ab 



This implies that II ([$]^) is symmetric for G L 1 ^. Moreover, if U in ( flOi) is symmetric, then 
the last identity in ([7]) holds. Again II is continuous, and II -1 real analytic, so that the proof of 
(ii) is completed. 

For case (iii), the third identity of ([5]) implies that the supplementary symmetry in ([8]) gives 



CL 



a,b G GL(L,C) 



a a 



b*b, r(ar 1 )*/ = a6 



This implies that /Il([$] rsj ) is skew-symmetric for [$]^ G L 1 ^. Again, if IU in 



TOD is skew- 
□ 



Let and denote the real analytic manifolds on the r.h.s. of the 

L,£ is obvious 



symmetric, then the last identity in ([8]) holds, completing the proof. 
Proof of Theorem [3l 

two equations in Theorem [3J Let us first show L 1 ^ = L 1 ^. The inclusion C 
because for a real representative $ the two conditions in ([7]) coincide. Moreover, this inclusion 
is continuous. Due to Theorem H|n) it is sufficient to show that LT^ 1 : U (L) PI Sym(L, C) — > L*, 
namely that one can choose a real representative in (flO|) . Let us use the fact that every symmetric 
unitary U can be diagonalized by an orthogonal matrix M G O(L), namely U = M t DM where 

e t6 L /2\ k e calculated with 



D = diag(e 101 , . . . , e l ° L ) with d x G [0, 2vr). Now let D\ = diag(e* 
the first branch of the square root and choose S = diag(o"i, . . . , cr^) G O(L) with oi G { — 1, 1} 
such that the phase of aie %e ^ 2 is in [0, it). Then let us introduce the unitary V = M t D^S. One 
has U = VV l . Furthermore set a = Jfe^) and b = —Qm(V) and then n _1 (t/) = [<3>]^ with 
a 
b 



Indeed II 1 is the inverse of II: 



n([<&]. 



7T 



V 
V 



u 



This construction of II -1 was done with a bit more care than needed, but it allows to show 
directly that II -1 is locally real analytic. Let E i— > U(E) be a real analytic path of symmetric 
unitaries. Then analytic perturbation theory [Katl Theorem II. 1.10] shows that the diagonal- 
ization U(E) = M(E) t D(E)M(E) can be done with analytic M(E) and D{E). Furthermore 
E i — ► D(E)z S(E) G diag(M/7rZ, . . . ,R/7rZ) with S(E) defined as above is also analytic because 
6 G R/2nZ i— > | G M/vrZ is analytic. Thus V(E) = M(E) t D^(E)S{E) is analytic and therefore 
also II -1 . 
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The proof of = L^; is just an adapted version of the usual rewriting of symplectic structures 
(here the symmetry induced by I and X) in terms of quaternions. Let the basis of EI (as real 
vector space) be 1 and the imaginary units k satisfying Hamilton's equations i 2 = j 2 = k 2 = 
ijk = —1. Then C is identified with the span of 1 and %. Now let us introduce 

T = ( o "0 1 1 j°i ) e Ma t(£x2L,H), 

where all blocks are of size § x §■ One readily verifies the matrix identity 

T*h/T = J - jX . (11) 

Hence one obtains a map T : Mat(2L xI,C) -> Mat(L x L, H) by matrix multiplication with 
T (from the left) which induces a map on the (complex) Grassmannians of right equivalence 
classes. Thus due to (11 ip we have exhibited an analytic map 

2.3 Symplectic group and Lorentz group 

Let K be one of the fields R, C and EI and let L be even if K = EI. The symplectic group SP(2L, IK) 
is by definition the set of complex 2L x 2L matrices conserving the Lagrangian structure in 
((7|) and ([8]) respectively, e.g. 

SP(2L, R) = { T € Mat(2L, C) | T*JT = J , T f JT = J } . 

One verifies that T 6 SP(2L,IK) if and only if T* 6 SP(2L,K). All symplectic matrices have a 
unit determinant. Using the Jordan form, it can be proven that SP(2L, IK) is arc-wise connected. 
Theorem [3] implies respectively the identity and isomorphism (direct algebraic proofs can be 
written out as well) 

SP(2L,M) = {T e Mat(2L, M) | T* R JT = J} , 

and 

SP(2L,H) = {T e Mat(L,H) | T* m IT = I} . 
More explicit formulas are given in the next algebraic lemma. 

Lemma 1 The complex symplectic group is given by 

A*C = C*A , A*D - C*B = 1 , B*D = D*B j . 

In this representation, elements o/SP(2L,IR) and SP(2L,EI) are characterized by having respec- 
tively real and self- conjugate entries A, B,C,D. 

As already became apparent in the proof of Theorems [3] and HI it is convenient to use the 
Cayley transform. The generalized Lorentz groups are introduced by 

V(L,L,K) = CSP(2L,K)C* . 



SP(2L,C) = U £ B D j e Mat(2L,C) 



S 



From the identities (jSJ) one can read off alternative definitions, e.g. 

U(L,L,R) = {T e Mat(2L x 2L, C) | T*QT = Q, T l J T = j) . 
Let us provide again more explicit expressions. 
Lemma 2 One has 

A* A - C*C = 1 , D*D - B*B = 1 , A*B = C*D j 

AA* - BB* = 1 , DD* - CC* = 1 , AC* = BD* j . 

Furthermore, in that representation, A and D are invertible and \\A~ 1 B\\ < 1 and ||.D _1 C|| < 1. 
For T G U(L, L, R) one, moreover, has C = B and D = A. For T G U(L,L,H) t/ie entries 
satisfy C = FBI and D = PAL 

Proof. The first relations are equivalent to T*QT = Q, the second ones then follow from 
the fact that T* G U(L, L, C) for T G U(L, L, C). The fact that A is invertible follows from 
AA* > 1. Furthermore AA* - BB* = 1 implies that A^BiA^B)* = 1 - A" 1 ^ 1 )* < 1, so 
that ||A _1 i?|| < 1. The same argument applies to D and D~ X C . The last two statements can be 
checked by a short calculation. □ 

2.4 Upper half-planes and Cartan's classical domains 

The upper half-plane and unit disc are defined by 

V C L = {Z e Mat(L x L, C) | i(Z* - Z) > } , = {[/ G Mat(L x L, C) | [/*[/ < 1 } , 

where Y > means that K is positive definite. Furthermore let us introduce the following subsets 
(here L does not need to be even for K = H): 

U£ = u£ n Sym(L, C) , = V C L n Self (L, C) , 

and 

B« = ro c n Sym ( jL) c), = n Self (L, C) . 

Let us note that ID? = n Skew(L,C). The sets 0£, Df and 7Df are called Cartan's first, 
second and third classical domain [Hua] . Furthermore Uf and ©f are also called the Siegel upper 
half-plane and the Siegel disc [Siej . The Cayley transform maps (via Mobius transformation) the 
upper half-planes bijectively to the generalized unit discs, as shown next. 

Proposition 1 The formulas 

u = (z - i i){z + zi)- 1 , z = % (i + u)(i - uy 1 , 

establish an analytic diffeomorphism from Uf onto Of for K = C, R, H. 



(-(L.I.C) = { ( ^ ^ ) g Mat(2I.€) 



A 5 



G Mat(2L,C) 
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Proof, (c/. |Sie] ; reproduced for the convenience of the reader.) If v G ker(Z + zl), then 
iv = —Zv so that < (v\i(Z* — Z)\v) = —2 (v\v) which implies v = 0. Hence Z + il is invertible 
and the first formula is well-defined. Similarly one checks the invertibility of 1 — U. To verify 
that one is the inverse of the other is a matter of calculation. Moreover, both formulas preserve 
the symmetry and self-duality of the matrices involved. □ 

The boundary <9Ug of Ug is a stratified space given as the union of strata cfyU^, I — 1, . . . , L, 
where diV^ is the set of matrices Z for which %{Z* — Z) > is of rank L — I. The maximal 
boundary is djJJ^ are the selfadjoint matrices. Furthermore <9Uf = fl Sym(L, C) with strata 
di\}\ = PI Sym(L, C). Corresponding formulas hold for K = EL 

Similarly the boundary <9D^ is a stratified space with strata <9;0^, I — 1, . . . , L, of matrices U 
for which U*U < 1 and rank(l — U*U) = L — I. For K = M, H one defines in the same way the 
stratified boundaries <9B^ = Ll/l^O 1 ^. Of particular importance will be the maximal boundaries 
(L even for K = EI): 

d L B^ = U(L) , d L Bl = U(L) n Sym(L, C) , d L 3^ = U(L) n Self(L, C) . 

By TheoremlHthe maximal boundary c^Of is hence identified with the Lagrangian Grassmannian 
. Let us also note that the Cayley transformation of Proposition [T] has singularities on the 
boundaries and mixes the strata. In particular, c^Uf is not mapped to djJDi^. 



2.5 Mobius action 

The Mobius transformation (also called canonical transformation or fractional transformation) 
is defined by 

T-Z = (AZ+B)(CZ+D)- 1 , T = ^ G GL(2L,C), Z G Mat(LxL, C) , (12) 

whenever the appearing inverse exists. This action implements the matrix multiplication, namely 
tt(T([$]„)) = T • tt([$]„) , if G Gr and [T$]„ G G£ v . (13) 

Indeed, let $ = ^ ^ j . Then [T$]^ G G£ v implies that Ca + Db is invertible. As 6 is invertible 

(because [$]„ G G£ v ), it follows that Cab" l + D = Cn([®]^) + D is invertible so that the Mobius 
transformation T • 7r([$]^) is well-defined. The conditions in f|T3|) are automatically satisfied in 
the situation of the following proposition. The proof of item (ii) is contained in the proof of 
Theorem below; item (i) then follows directly from (ii) due to Proposition [TJ 

Proposition 2 [Hual ISTe] Let K = C, R, H and let L be even for K = H. 

(i) SP(2L, K) acfo on by Mobius transformation. 

(ii) U(L, L, K) acts on Of 6y Mobius transformation. 

The following proposition states that the action of Proposition [2]Jn) extends to the stratified 
boundary of D^. Moreover, the action on the maximal boundary c^Df implements the natural 
action of the symplectic group on the Lagrangian Grassmannian. Due to singularities it is not 
possible to extend the action of Proposition ^i) to any stratum of the boundary of Uf . 
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Theorem 5 Let IK = C, R, H and L even if K = H and I = 1,...,L. The Lorentz group 
U(L, L, IK) acts on 80^ by Mobius transformation. For the case I = L of the maximal boundary, 
this action implements the action o/SP(2L, K) on the Lagrangian Grassmannian L^: 

n([T$].) = crc*-n([$].), Tg SP(2L,K). 

Proof. One has to show that for U G c^HDf and T, T' G U(L, L, C) the Mobius transformation 
T • U is well-defined, is again in <9;Of and that (TT f ) ■ U = T ■ (T 7 ■ [/). Let T be given in 
terms of A, B, C, D as in LemmaEl Then this lemma implies that (CU + D) = D(l + D~ 1 CU) is 
invertible so that the Mobius transform is well-defined. Let us first show that (T -U)*(T -U) < 1. 
For this purpose, one can appeal to the identity 

(C U + D)* (C U + D) - (AU + BY (AU + B) = 1 - U* U , (14) 

following directly from the identities in Lemma [2j Indeed, multiplying fTl4|) from the left by 
(CU + D)*)- 1 and the right by (CU + D)- 1 and using 1 - U*U > shows (T ■ U)*(T ■ U) < 1. 
Next let us show that the invertible (C77 + D) maps ker(l - U*U) to ker(l - (T ■ U)*(T ■ U)) 
and ker(l - U*Uy to ker(l - (T ■ U)*(T ■ U)) x . If v Gker(l - U*U), then (JUJ implies 

|| (CU + D)v\\ = \\(AU + B)v\\ = \\T-U(CU + D)v\\, 

so that (CU + D)v G ker(l - (T ■ U)*(T ■ U)) because 1 — (T ■ C/)*(T • C/) > 0. Similarly for u G 
ker(l-[/*[/) ± onehas ||(CC/+D)u|| > ||T-f/(Cf/+ J D)t;|| implying that v £ ker(l- (T-U)*(T-U)). 

The argument up to now shows that T ■ U G cPJ^. A short algebraic calculation also shows 
that (TT 1 ) ■ Z = T ■ (T' ■ Z). It remains to show that the symmetries are conserved in the cases 
K = R, EL For T G U(L, L, R) one has C = 5 and D = A by Lemma El so that 

T-U-(T-UY = ((BU + A)- 1 ) 1 [(BU + 1)\AU + B) - (AU + B) l (BU + A)] (BU + A)' 1 . 

For symmetric U one checks that the term in the bracket vanishes due to the identities in 
Lemma El implying that T ■ U is again symmetric. Similarly one proceeds in the case K = EL 

Now let us come to the last point of the proposition. Given U G di0f;, let [$]^ G be 
such that Il([$]^) = U (by the construction in the proof of Theorem HJ). For T G SP(2L,K), 
one then has [T$]„ G Lf. Theorem H implies that both [C$]~ and [CT$]^ = [CTC*C$]^ are in 
Gjf. From flTSJ) now follows that CTC* ■ 7r([C$]^) = ?r([CT$]^). □ 

It is interesting to note (and relevant for Section 13.61) that the Mobius transformation sends 
to Uf" for some matrices in GL(2L,C) which are not in SP(2L,K). In particular, for 5 > 
one even has: 

( J ^ Yz = Z + i5l G Uf, for ZGU^U^. (15) 

The following formula shows that the Mobius transformation appears naturally in the cal- 
culation of a volume distortion by an invertible matrix, namely the so-called Radon-Nykodym 
cocycle. It will be used for the calculation of the sum of Lyapunov exponents in Section 14.11 
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Lemma 3 Suppose that [$]^ G G£ v and [T$]^ G G£ v where T G GL(2L, C). With notations 
for T as in (1120 one /ios 



det((T$)*(T$)) det((r-7r([^))*(T-7r([$]^)) + l) 
det($*$) ~ det((7r([$]„))*(7r([$]„)) + l) 



det(C +D) 



a 



Proof. Let $ = ( , ] . As b is invertible by hypothesis, 



det( (T$)*(T$) ) det( {Aab- 1 + B)* {Aab~ l + B) + (Cab~ l + £>)* (Cao" 1 + D) ) 
det($*$) det( (afe- 1 )* (afe- 1 ) + 1) 

Furthermore, it was supposed that CoJcr 1 + D is also invertible, and this allows to conclude the 
proof because 7r([$]^) = ab^ 1 . □ 

2.6 Singular cylces 

Given £ G Sym(L, R) n Self(L, C), let us set = ^ ^ V The associated singular cycle (often 
called Maslov cycle, but it actually already appears in Bott's work [Botj) is 

= G $ c L n ^ C C L ^ {0}} . 

It can be decomposed into a disjoint union of h L , / = 1, . . . , L, where 

L^' z = {[$]„ G | dim($ C L n ^ ? C L ) = /} . 

It is possible to define singular cycles associated to Lagrangian planes which are not of the form 
of but this will not be used here. 

It is convenient to express the intersection conditions in terms of the Wronskian associated 
to two 2L x L matrices $ and \I> representing two Lagrangian planes 

namely one checks that $ C L n ^ C L ^ {0} <=> det(W = 0, and, more precisely, 

dim($C L n ^C L ) = dim(ker^($,^)) . (16) 

Note that even though VT(<&, VP) does depend on the choice of basis in the two Lagrangian planes, 
the dimension of its kernel is independent of this choice. Furthermore, W(<&, = if and only 
if = 

Arnold showed in |Arnlj that L^'^ is two-sided by exhibiting a non- vanishing transversal vector 
field on . This allows to define a weighted intersection number for paths in a generic position, 
namely for paths having only intersections with the highest stratum L^'^' 1 . Bott proved a similar 
result for L 1 ^ already earlier |Botj . Using the following proposition, it will be straightforward in 
the next section to define intersection numbers for paths which are not necessarily in a generic 
position. 
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Proposition 3 11(17^'^) = d^ l 3f where c^Df is the following subset of the maximal boundary 
L Bf: 

d^B^ = {Ued L B^ | rank(II([^]„)-l7) = L - 1} . 

Setting d|Bf = U l=lt ... :L d^ l Bf, one has n(L^) = dfDf. One /las 1 £ dfOf / or an ?/ £ e 
Sym(L,R)nSelf(L,Kj.' 

Proof. Given [/, let $ = ^ ^ J where a = \{U + 1) and 6 = |(C/ — 1) as in the proof of 
Theorem |H Hence II ([$]^) = £7. Then one verifies 

w{*t,s) = a-tb = i + i -,£(17-1)) = i(i-«e)(r/-n([^)). 

As 1 — is invertible, the result follows directly from (jT6l) . No more care is needed in the real 
case because the dimension of the kernel of the Wronskian is independent of a basis change from 
the above $ to the real one in the proof of Theorem HI □ 

We will only use the singular cycles associated to £ = — cot(^) 1 where ip G (0, 2tt), and 
only need to use K = C as the supplementary symmetries are irrelevant for the definition and 
calculation of the intersection number (Bott or Maslov index) in the next sections. Let us denote 
these singular cycles, subsets of 17^, by Vf L = U^i^.^L^'. The image under II will then be 
written as <9£0^ = U^i^.^c^O 1 ^. Because 1I(\I/^) = e lip l, it follows from Proposition [3] that 

n(L^) = = { U G U(L) | e'f eigenvalue of U with multiplicity 1} . (17) 

2.7 The intersection number (Bott index) 

Let T = {[§ E ]rJ) Ee[E ,Ei) be a (continuous) closed path in I7£ for which the number of intersections 
{E G \Eq,Ei) I T(E) G I7£} is finite. At an intersection T(E) G let O^E'), . . . , 6 t (E') be 
those eigenphases of the unitary U(T(E')) which are all equal to ip at E' — E. Choose e, 5 > 
such that 9k(E') G [ip — 5, ip + 5} for k = 1, . . . , I and E' G [E — e, E + e] and that there are no 
other eigenphases in [ip — 5, ip + 5} for E' ^ E and finally Ok(E') ^ ip for those parameters. Let 
n_ and n + be the number of those of the I eigenphases less than ip respectively before and after 
the intersection, and similarly let p_ and p + be the number of eigenphases larger than ip before 
and after the intersection. Then the signature of T(E) is defined by 

sgn(r(77)) = i (p + -n + -_p_+n_) = Z-n+-p_. (18) 

Note that — / < sgn(r(77)) < / and that sgn(r(77)) is the effective number of eigenphases that 
have crossed ip in the counter-clock sense. Furthermore the signature is stable under perturbations 
of the path in the following sense: if an intersection by L^'' is resolved by a perturbation into 
a series of intersections by lower strata, then the sum of their signatures is equal to sgn(r(77)). 
Finally let us remark that, if the phases are differentiable and dE@k{E) ^ for k — 1, . . . , I, then 
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sgn(r(£J)) is equal to the sum of the I signs sgn(dE0k(E)) , k = 1, . . . ,/. Now the intersection 
number or index of the path T w.r.t. the singular cycle Vf L is defined by 

ind(r,LD = ^n(T(E)). (19) 

r(E)ehl 

Let us give a different expression for this index. If E G [E , E{\ \— > 9f are continuous paths 
of the eigenphases of U(T(E)) with arbitrary choice of enumeration at level crossings, then each 
of them leads to a winding number. A bit of thought shows that 

L 

ind(r,L£) = Wind( E G [£ ,£i) 0? ) . (20) 
i=i 

In particular, the r.h.s. is independent of the choice of enumeration at level crossings. Moreover, 
for the r.h.s. to make sense, one does not need to impose that the number of intersections 
is finite, as is, of course, necessary in order to define an intersection number. Similarly, the 
index of a closed path T in the real or quaternion Lagrangian Grassmannian could be defined; 
however, this index coincides with ind(r, L£) if T is considered as path in the complex Lagrangian 
Grassmannian. 

In the literature, ind(r,L^) is often referred to as the Maslov index, at least in the case of 
a path in L*. The same object already appears in the work of Bott [Bot] though, and it seems 
more appropriate to associate his name to it. The above definition using ffT8]) appears to be 
considerably more simple, and the author does not know whether it was already used elsewhere. 

2.8 Arnold's cocycle 

Arnold has shown [ArnT] that # x (Lf , Z) = Z. This and H\I% , Z) = Z follows from Theorem H 
The generator uj of the de Rahm groups can be chosen as follows. A continuous closed path 
T = ([^j^EelBo.-Ei) in gives rise to a continuous closed path E G [E ,Ei) i-> det(n([$ £; ]^)) 
in S 1 . Its winding number defines the pairing of (the de Rahm class of) uj with the (homotopy 
class of the) path Y: 

(uj\T) = Wind( E G [£ ,£i) det(n([$ £; ]^)) ) . 

Any continuous path in hj] can be approximated by a different iable one, hence we suppose from 
now on that T is differentiable. Then one can calculate the pairing by 

f El dE 

(u\T) = —%md E log(det(n([$ s ].))) . (21) 
Je 2tt 

The next theorem is Arnold's main result concerning this cocycle. 

Theorem 6 |Arnl] Provided a closed differentiable path T has only finitely many intersections 
with the singular cycle Vf L , 

ind(r, hi) = ( w | r ) . 
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Proof. Set U(E) = U(T(E)). As already used in the proof of Theorem HJ the diagonalization 
U(E) = M(E)*D(E)M(E) can be done with a differentiable unitary matrix M(E) and a differ- 
entiate diagonal matrix D(E) = diag(e<, . . . ,e<). As M E (d E M E )* = -(d E M E )(M E )\ one 
has 

i L 
3m d E log(det(tf(£))) = - Tr ((D(E))*d E D(E)) = V d E 9 E . 

Integrating w.r.t. E as in (121 p hence shows that (u;|r) is equal to the sum of the winding 
numbers of the eigenphases and this is equal to the index by fl20l) . □ 

Remark It follows from the Gohberg-Krein index theorem that the intersection number is also 
equal to the Fredholm index of an associated Toeplitz operator [BSJ. 

If T G SP(2L, C), then one can define another closed path in by TV = ([T§ E ]^) Ee [ EojEl y 
Because SP(2L,C) is arc-wise connected, it follows from the homotopy invariance of the pair- 
ing that ( lu | Tr ) = ( u) | r ) . The following proposition allows T to vary and also analyzes 
intermediate values of the integral in fl2T|) , denoted by: 



Id 



E 



de 
2^ 



3m d e log(det(n([$ e ]^))) , E G [E , E 1 



Proposition 4 Let T = ([<& E ]~) e^EcEx) be a closed differentiable path in and (T E ) E( z[ Eo , El ) 
be a differentiable path in SP(2L, C) such that V = ([T E &]^) Ee[E ,Ei) is a closed path in for 
any given [&\~ G h L . Furthermore let us introduce the closed path T" = ([T e <& e ]^)e£[Eo,Ei)- 
Then 



Hr") = mo + (u\V) . 



Furthermore, with the notation 



cr E c* 



A E B E 

C E D E 



one has {independent of [$]~) 
and, uniformly in E, 

-E 



Eo 



dE 
2^ 



3m d E log(det(A B ( J D E )- 1 )) 



(22) 



< L 



Proof. Set U e 

-E 



CO 



f a de 

oj - —% m d e Wdet^D 6 )- 1 )) 
Je 2vr 

Due to Theorem [5], 

/ _£ Q m 8 e ( log(det(A e U e + B e ))- log(det(C e U e + D e )) 
J E n 2tt V 



n([$ 

-E 
'E t) 
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In the first contribution, let us use det(A e U e + B e ) = det(B e (U e )* + A e ) det{U e ), of which the 
second factor leads to jjf uj. Hence 



f E r E t m de 

/ lo - to — I — 3m d e (log(det(B e (f/ e )* + A e )) - log(det(C e f/ e + D e ))) . 
Jr" Jr Jeo 2tt 



>E 

Now A e and D e are invertible due to Lemma [21 Hence 

E pE pE 

UJ — I UJ 

v" Jr je c 



de 

/ — 3m «9 e log(det(A e ( J D e )- 1 )) 
Je 

/ J- Q e 3m ( Tr (log(l + (A e )~ 1 B e (U e )*)) - Tr (log(l + (D e ) _1 C e U e )) ) 

J Ea 27T V / 



'E 

As || (A e ) _1 i? e || < 1 and || (D e ) _1 C e || < 1 by Lemma [21 only one branch of the logarithm is 
needed. As E — > E\ this term therefore vanishes implying the result on the winding numbers 
because the third term on the l.h.s. then gives ( u \ V) as one sees repeating the above arguments 
with U e replaced by the constant n([$]^). Moreover, one may carry out the integral on the r.h.s. 
of the last equation using the fundamental theorem, so that this r.h.s. is equal to 

^ ( Tr (log(l + (A E )' 1 B E {U E T)) - Tr (log(l + {D E y l C E U E )) ) . 

The bound follows now from the spectral mapping theorem for the logarithm function. □ 
In order to calculate the integral in fT2T]) . one can also appeal to the following formula. 

Lemma 4 If E i— > & E is differentiable and [& E ]~ G then 

3m 8 E log(det(n([$ B ]^))) = 2Tr( J (d E . (23) 

Proof. As in the proof of Theorem El let us begin with the identity 

id E log(det(H([$^))) = Tr( (HQ^))* (d E U([^)) ) . (24) 

Introducing the invertible L x L matrices 

<P E = (1 ±zl)$ E , 

one has nQ^]^) = 0^(0^) _1 . Furthermore the following identities can be checked using the 
fact that <& E is Lagrangian: 

(</>+)>+ = {4> E -)*4>- = ($ B )*<^, {^lTd E <j) E = (<S> E )*d E $ E ± i($ E )*Jd E <S> E . (25) 
Corresponding to the two terms in 
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there are two contributions G\ and Ci in ([21]). For the first one, it follows from cyclicity 

while for the second 

C2 = -Tr(^)- 1 ^^) = -iv(((^)^)- 1 (^r^0j). 

Combining them and appealing to (1231) concludes the proof. □ 

Next let us calculate the pairing for two examples. The first one was already given by Arnold 
[Arnlj . For 77 G [0, 7r], introduce the symplectic matrices 

cos(?7) 1 sin(?7) 1 \ / e"? 1 



^ V-sin(r7)l cos(r/)l J ' C Rr,C * ' \ 1 

as well as the closed path Y = (i2 J7 [$]^) 7?6 [ 0i7r ) for an arbitrary [$]^ G L£. As det(n([_R^$]^)) = 
e 2iLv det(n([$]^)), one deduces (u\T) = L. The second example concerns transfer matrices. 

Lemma 5 Let T E be a matrix built as in (|2J) from a selfadjoint matrix V and positive matrix 
T. For an arbitrary [$]^ G with (10)$ invertible and M. GSP(2L,C), consider the path 

r = (mt e [^) 

E£R where K — RU{oo} is the one-point compactification. Then T is closed and 
(oj |T) = L. Moreover, E G R 1— » oj is strictly monotonously increasing. 

Proof. First suppose that M = 1. Then note that Il([T jE $]^) = 1 + 0(E' 1 ), hence the path 
is closed. Next one calculates 

E If (El - V)T- 1 -i{T + T- 1 ) (El - V)T- l +t(T-T- 1 ] 

2 V (El - V)T~ X -i(T-T- 1 ) (El - V)T~ l + z(T + T~ 1 ; 

Let A E and _D B denote the upper left and lower right entry. Then d E log(det(v4' E )) is equal 
to 

Tr ((T-^E 1 - VfT- 1 + (T + T" 1 ) 2 + i(T~ l VT - TVT- 1 ))' 1 ^ + T~ 

As there is a similar expression for Sm d E log(det(D £; )), this shows that the integral in (122]) is 
finite. In order to calculate the winding number, it is convenient to consider the homotopy 

- ( iE1 "n^ r(A) "' "?')■ 0£A£1 - 

where V^(A) = AV^ and T(\) = AT + (1 — A)l (the latter is always positive). Then the pairing of 
T(A) = (T E (\)[§}^) Ee n with oj is independent of A. Hence it is sufficient to calculate the pairing 
at A = 0, which is due to the above replaced into ( 122]) given by 



/CO JP 
— Tr(( J B 2 + 4)- 1 2 1) = L 
-oo 
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completing the proof in the case M. = 1. If M ^ 1, the winding number is the same because 
SP(2L, C) is arc- wise connected, as already pointed out above. The monotonicity can be checked 
using the r.h.s. of ([23]) for $ £ = MT E §. In fact, is strictly positive, and 

($ E yjd E $ E = <r ( T Q 2 ° ) $ , 

which is strictly positive by hypothesis. As the trace of a product of two positive operators is 
still positive, this concludes the proof. □ 



3 Jacobi matrices with matrix entries 

Given integers L, N G N, let (T n ) n= 2,...,jv and (V n ) n= i^„ t N be sequences of respectively positive 
and selfadjoint Lx L matrices with complex entries. Furthermore let the left and right boundary 
conditions ( and £ be also self-adjoint L x L matrices. In the real and quaternion case, one 
chooses £ and £ symmetric and self-dual. Then the associated Jacobi matrix with matrix entries 
is by definition the symmetric operator acting on states (ft = (<p n )n=i n £ ^ 2 (1, • • • , N) (g> C L 

by 

{Hf4>) n = T n+1 (j) n+l + V n (f> n + T n <\) n _ x , n = 1, . . . , N , (26) 
where T\ = TV+i = 1, together with the boundary conditions 

00 = C 01 ) 4>N+1 = £ 4>N ■ (27) 

One can also rewrite defined in (I26p and (j27j) more explicitly as a block diagonal matrix; 
this gives ([T]) albeit with V\ and replaced by V\ — ( and Vjy — £• The dependence on £ is not 
specified, but it could potentially be used for averaging purposes. If ( = £ = one speaks of 
Dirichlet boundary conditions. 



3.1 Transfer matrices and dynamics of Lagrangian planes 

As for a one-dimensional Jacobi matrix, it is useful to rewrite the Schrodinger equation 

Hg<f> = E<j>, (28) 

for a complex energy E in terms of the 2L x 2L transfer matrices T E defined in (j2J). For a real 
energy E el, each transfer matrix is in the symplectic group SP(2L, C). If H is, moreover, real 
or self-dual, then the transfer matrices are in the subgroups SP(2L, M.) and SP(2L, H) respectively. 
The Schrodinger equation fT28|) is satisfied if and only if 

and the boundary conditions (|27|) hold, namely 

f T ^ ) G $ f C L , f TiV+ ^ iV+1 ) G , (29) 
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where we introduced for later convenience the notations 




Both of the two L- dimensional subspaces of C 2L appearing in the conditions (129!) are Lagrangian. 
One way to search for eigenvalues is to consider the whole subspace in the left equation of ( 1291 ). 
then to follow its evolution under application of the transfer matrices, and finally to check whether 
at N the resulting subspace has a non-trivial intersection with the subspace on the r.h.s. of 
( 1291) . For perturbation theory in Section [U it is useful to incorporate a symplectic basis change 
M. G SP(2L, C) which can be conveniently chosen later on. In the one-dimensional situation 
this corresponds to pass to modified Priifer variables. If H N is real or self-dual, one chooses 
Ai eSP(2L,IR) or M. G SP(2L, H). As above, half-dimensional subspaces will be described by 
2L x L matrices <3>^ of rank L composed of column vectors spanning it. Then their dynamics 
under application of the .M-transformed transfer matrices is 

$f = MT^M-^ti, *g = M* c . (30) 

If E G M., these planes are Lagrangian. As the boundary condition on the left boundary is satisfied 
automatically (it is chosen as the initial condition), the second condition in (12"§|) multiplied by 
M. together with the Wronski test ([TBI) leads to 

multiplicity of E as eig envalue of = L - rank ( (M ^ 5 )* J . (31) 

This implies also 

Proposition 5 Let T = ([$f]~) SeS . For every £ G Sym(L, R) n Self(L,C), the set {E G 
K. | T{E) G L^'^} of intersections with the singular cycle L^'^ is finite. 

The dynamics (1301) is more easily controlled under the stereographic projection. Let us first 
consider the case ^sm(E) > 0. In this situation the stereographic projection tt of (!30j) gives a 
dynamics in the upper half-plane U^, or Uf if H is real and if H is self-dual. In fact, 
Zf = 7r($f ) = A4 ■ (E 1 — V\ — Q is in U 1 ^ and, moreover, the transfer matrices factor as follows: 

E (1 i$m(E)l\ ( ($te(E) 1 - V n ) T* 1 -T n \ 

n V 1 A ^ J ' 

The matrix on the right is in the symplectic group SP(2L, C) acting on U^, the one on the left 
also sends U L to U L because of 05). The same applies for real and self-dual H N . Hence the 
following Mobius action is well-defined: 

Z E n = MXfM-'-Z^, Zf = M-{E1 -V x - C), (32) 

In the case M. = 1, this is a matricial Ricatti equation Z E = E — V n -\-T n {Z^_^)~ l T n . Comparing 
with PJD, 

Z° = n = l,...,N. 
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In particular, [$^]~ is in the domain G£ v of 7r whenever Qm(E) > 0. If the boundary condition 
( is invertible, one may also set Z E = A4 ■ and then Zf = MT E M.~ l ■ Z E , even though Z E 
is not in U^J. Furthermore, the map E G TJf ^ Z§ G is analytic (Herglotz). The map has 
poles on the real axis as can be read off the Dean-Martin identity (T^oT) proven below, but will 
not be used in the sequel. 

As real energies are not always permitted, the Z E are not convenient for the calculation of 
the eigenvalues. For any E with ^sm(E) > 0, let us rather use 

u E = n([<&fu, n = 0,...,N. 

For real E, this is well-defined and U E is a unitary because of Theorem HI This unitary is 
symmetric or self-dual if H is real or self-dual. Iterating (I30I) and recalling the definition of 
the stereographic projection shows that U E is actually of the explicit form given in ([3]) if one 
chooses M. — 1. Hence this proves Theorem [T^i) and part of (ii). For Qm(E) > 0, the above 
arguments imply Z E is well-defined and hence also U E = C ■ Z E . One, moreover, concludes that 
U E = C ■ Z E is in the generalized unit disc (for n ^ 0). The dynamics is given by 

U E = CMT E M~ X C -XJ E _^ U E ^•(l-^CMl + ^r 1 - (33) 

For Qm(E) > 0, this is just the Cayley transform of ( 1321) . while for E G R, it is the dynamics of 
Theorem [5j The following lemma proves Theorem [H(ii) and the first part of (iii). 

Lemma 6 The map E i— ► U E is analytic in a neighborhood of U£ = Uf U dUf. At level 
crossings, the eigenvalues and eigenvectors can be enumerated such that they are analytic in a 
neighborhood o/Uf as well. 

Proof. Analyticity of U E away from the real axis follows from the analyticity of Z E because 
U E = C ■ Z E and the inverse in the Mobius transformation is also well-defined, cf. Section 12.41 
Moreover, the characteristic polynomial is a Weierstrass polynomial that has a global Puiseux 
expansion which is analytic in the Lth root of E. Hence the eigenvalues and eigenvectors can be 
chosen (at level crossings) such that they are analytic in the Lth root of E (e.g. |Katl Chapter 
II]). It will follow from the arguments below that the Puiseux expansion actually reduces to a 
power series expansion in E. 

Now we analyze in more detail the situation in a neighborhood of the real axis. The plane 
$fr is a polynomial in E. Let us use the notations $^ = ( ^ E J. It follows from the argument 

in (Q that a E + ih E has maximal rank for E G R so that £g!h det(a B + ib E ) has no zero. 
Moreover, one has the large E asymptotics 

Qfi = E N M ^ n^i^n 1 ^ + 0{E N - 1 ) . (34) 

which implies 

det(a E + ib E ) = E NL det [( 1 j V M (l )) JJdet^)- 1 + 0(E NL ~ 1 ) . 



n=l 
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Therefore infe G R | det(a E + ib E )\ > and the infimum is realized at some finite E. A perturbative 
argument shows that also inf Ees s I det(a E +i b E ) \ > where Sg = {E G C | \$sm(E) \ < 5} is a strip 
of some width 5 > 0. Calculating the inverse with the Laplace formula shows that (a E + tb E )~ 1 
is also analytic in S$- Thus also U E = (a E — i b E )(a E + % b E )~ 1 is analytic in S$ and thus analytic 
in a neighborhood of due to the above. 

Finally one can appeal to degenerate perturbation theory [Kat[ Theorem II. 1.10] in order to 
deduce that the eigenvalues and eigenvectors of the unitary matrix U E (hence E real) are also 
analytic in a neighborhood of the real axis, that is are given by an analytic Puiseux expansion. 
As this neighborhood has an open intersection with the upper half-plane, the above Puiseux 
expansion is therefore also analytic, namely only contains powers of E. □ 

It follows from (JMD that 

U E = CMC* 1 + 0{E- 1 ) . (35) 

Let < di < 2n be the eigenphases of the symmetric unitary CM.C* ■ 1. The eigenvalues of 
U§, chosen to be real analytic in E G R as in Lemma El are denoted by e i6,AU , I = 1, . . . , L. The 
eigenphases are chosen such that 8 E l — > Of 4 for E — > — oo. In the case M. — 1, one hence has 
6 E t — > for E — > — oo as in Theorem [lj 

Let us conclude this section by choosing particular right boundary conditions, namely, for 
^G(0,2vr), 

£ = - cot(|) AT 1 ■ 1 =► [M9t]„ = [*-a*mi]„- 

The corresponding Hamiltonian will be denoted by H^. If M. = 1 and <p — 7r, these are Dirichlet 
boundary conditions on the right boundary. Due to Proposition [3] and ffTTl) . the eigenvalue 
condition fl3Tl) becomes 

multiplicity of E as eigenvalue of = multiplicity of e %v as eigenvalue of U§ . (36) 

Setting Ai = 1 and if = tt, this proves Theorem QJry). 

3.2 Monotonicity and transversality 

This section provides the proof of Theorem [Tfv) (just set M. = 1 in the below). Of course, 
the second statement of Theorem D^v) follows immediately from the first one upon evaluation in 
the eigenspaces of U E . The following proposition also shows that the curve T = #el is 

transversal to the singular cycle and always crosses it from the negative to the positive side. 

Proposition 6 For E G M and N > 2, one has 

- i (u E )*d E u E > 0. 

Proof. As in the proof of Lemma HI let us introduce </>± = ( 1 ±*1 ) Qff- These are invertible 
Lx L matrices and one has U§ = <p^{<f)f ) _1 = ((^) _1 )*(0f )*. Now 

(u E yd E u E = 



)*d, 



)*d f 
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Thus it is sufficient to verify positive defmiteness of 

- \{&d E <p E - {<j>%Yd E <i>i\ = 2 {$*yjd E *% . 

From the product rule follows that 

N / N \ /n-1 \ 

a B *% = y: m n ^ (^) n ^ $ c 

n=l \Z=n+l / \/=l / 

Using M*JM = J, this implies that 

N /n-1 \ * /n-1 

M)*^&<^ = e $ c n ^ ( t « £ )^ ( d * r n) n ^ ) $ c • (37) 

n=l \i=l / \l=l 

As one checks that 



(T^)- 1 


each of the summands in (1371) is positive semi-definite. In order to prove the strict inequality, it 
is sufficient that the first two terms n — 1, 2 in fl37j) give a strictly positive contribution. Hence 
let us verify that 

CT(<rs)? + ( (T rs)>- 

For this purpose let us show that the kernel of the matrix on the l.h.s. is empty. As ({T 2 E )*) 1 = 
—JT E J, we thus have to show that a vector ( ^ j e C 2L satisfying 



(TO*)- 1 o \ E ( v \ s / (TO)- 1 o W v 

0i i2 l«;j~ i2J l U 



actually vanishes. Carrying out the matrix multiplications, one readily checks that this is the 
case. □ 

Of course, one can regroup the terms in (1371) into packages of two successive contributions and 
each of them is positive by the same argument. If this bound is uniform for the packages (e.g. 
the spectrum of the T n is uniformly bounded away from 0), one actually deduces an improved 
lower bound by C E N for some C E > 0. 

3.3 The total rotation number 

In this section, we complete the proof of Theorem (TJ in particular the second part of item (iii). 
Throughout Bel. The total rotation number is defined by 



/a 
de 9m d e log(det(L/^))) . (3* 
-oo 
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It will be shown shortly that the integral converges. Using the notations of Section I2.8[ Gj| = 
Jjf uj for T = ([^jvWsel- Let 6f il be the analytic eigenphases of U§ as introduced after (1351) . 
We deduce after integration of ( 1381) that 

= E Ki - en • (39) 

1=1 

This justifies the term total rotation number. The following result could also be deduced from 
the results of the previous section, but its proof (a homotopy argument) directly completes the 
proof of Theorem [T^iii). 

Proposition 7 The total rotation number 0j| is well-defined and satisfies 

lim 6f = 0, lim 6^ = 2n N L . 

E— >— oo i?— >oo 

Proof. Expanding $f = M ]!n=i $ C shows 

((*5)**5) _1 = U E2N ( J o ) N M * M ( J o + 0{E2N ~" 



E~ 2N ( J J M*M( I ) + 0(£T 2 - V - J ' 



Similarly one verifies 
Hence follows 



so that the integral in fl38l) exists due to ff23l . which can alternatively be derived from Lemma [51 
Furthermore, one deduces from fl35|) that T = ([^tvWegS is a closed path in and its winding 
number is given by 

Mr> = -L lim e*. 

In order to calculate the winding number and prove the last statement of the proposition, one 
applies Proposition H] and Lemma [5] iteratively to the path T = {[M. Y\n=i Xf^cl^EeW- Alter- 
natively one can use a homotopy H N (X) from H N (1) = to H N (0) which is the sum of L 
un-coupled one-dimensional discrete Laplacians (using the homotopy of the proof of Lemma [5] on 
every site n). For each of the one-dimensional discrete Laplacians the winding number is again 
easy to calculate and equal to N. □ 

Proof of the last statement of Theorem [T](iii) . The homotopy discussed at the end of the proof 
of Proposition [7J is analytic and hence one deduces 9fj l — — > 2irN for E — > oo for each / as 
this is the case for the one-dimensional discrete Laplacian. □ 

Remark Using Proposition [7J one can also give a nice alternative proof of 8e0§i > 0. Indeed, 
according to fl3T)|) . the unitary U§ can be used in order to calculate the spectrum of H£ for every 
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if G (0, 2tt). Counting multiplicities, this spectrum consists of NL eigenvalues. By Proposition [7] 
and Theorem [6j the total number of passages of eigenvalues 8§ L by ip (intersections with the 
singular cycle L£) is bounded below by NL. As there cannot be more than NL, all these 
passages have to be in the positive sense because a passage in the negative sense would lead to 
at least two more eigenvalues of H N . 



3.4 Telescoping the total rotation number 

By the results of the last section and due to the fact that a change of right boundary condition 
can change the number of eigenvalues by at most L, one has 



2tt 



9^ - # { eigenvalues of iff < E } 



< 2L 



(40) 



Hence 0^ allows to count the eigenvalues of H N up to boundary terms. For this purpose it is 
useful to telescope Q E into N contributions stemming from the L- dimensional slices: 



u v 



N 

E 

n=l 



de 3m d e log 



detgg 
det(^_^ 



Here we have used ( 1331) and the fact that Uq is independent of e. This is indeed a good way to 
telescope because = C M. T£ M.~ l C* ■ U^_ x so that Proposition @] and Lemma [5] imply that 
each summand satisfies 

-E 

de 3m d e log 



( &et{C MT£ M- 1 C* -V e n _Z 



det(^_ i ; 



< 2L 



(41) 



Moreover, with the notation 

CMT r fM- l C* = 
the same calculation as in Proposition H] implies that 

N „E 

©f = E / de ^ m d * lo g( det ((^. + B* n {UUT){D e n + ClUUY 1 )) ■ (42) 

i J — oo 



n=l 



Now it is actually possible to apply the fundamental theorem in every summand by determining 
the branch of the logarithm uniquely from the transfer matrix T E , and independent of the prior 
transfer matrices. Indeed, one can factor out det(A^(D^) _1 ) and use the fact that (A^j^B^ and 
{D^)~ l C^ have norm less than 1 by Lemma [2J so that as in the proof of Proposition HJ 



U 7V 



JV 

E 

n=l 



3m 



de d e log(det(^( J D:)- 1 )) 



(43) 



+ Tr (log (1 + {A E )~ X B E (U^T )) - Tr (log (l + (D E )^ C E U E _, )) 
A refined version of this formula is exploited in Section |H 
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3.5 The Dean-Martin identity 

In this and the next two sections, M. = 1 and £ = (Dirichlet boundary conditions on the right 
boundary). The Hamiltonian will be denoted by H N without a further index. According to (I3TT) . 
E is an eigenvalue of H N if and only if 

det ((0 1)J$%) = det (V 0) l[T n E $ c J = . 

Due to lj3lj) . the multiplicity of the zero is the multiplicity of the eigenvalue. Therefore the l.h.s. 
of this equation is a polynomial of degree NL in E with zeros exactly at the NL eigenvalues 
of H N . Comparing the leading order coefficient, one deduces a formula for the characteristic 
polynomial: 

TV 

det{El-H N ) = det ((l 0) $f ) JJdet^ 1 ) . (44) 

71=1 

In order to find a recurrence relation for the characteristic polynomials, let us suppose that 
%m(E) > and note that Z§ = (l 0) $f ((l 0) ^n-i)~ 1t n- Taking the determinant of this 
formula, the identity (j4"4l) applied twice gives 

det (Z§) = ^ -4r . 45 

v w; det (El -H N - 1 ) v ; 

Let us call this the Dean-Martin identity, due to the contribution [DM] . These authors then used 
the identity (I4"5"j) at real energies in order to calculate the spectrum of H N_1 by counting the 
singularities of det(Z^). This can be made more explicit by adding a small imaginary part 5 > 
to the energy. Then consider the path E G M i— > det(Z E+tS ) G C. Even though Z E+tS is in the 
upper half-plane, its determinant may well have a negative imaginary part. However, it never 
takes the values ± i. Now each passage of the path (near) by an eigenvalue of H 1 ^ -1 leads to an 
arc in either the upper or lower half-plane with passage by either i/5 or — 1/8, pending on the 
sign of the numerator in (j4"5l) . Both arcs turn out to be in the positive orientation. A multiple 
eigenvalue leads to a multiple arc. The topologically interesting quantity is the winding numbers 
of the path around i and —i. Calculating the sum of the corresponding phase integrals gives a 
total rotation number which in the limit 5 —>■ coincides with 0^. This allows to give a nice 
alternative, but considerably more complicated proof of Proposition [7J In the case L — 1, all the 
arcs are in the upper half-plane and the argument just sketched is particularly simple. 

3.6 Green's function and continued fraction expansion 

The aim of this short section is to illustrate the use of the dynamics in the upper half-plane. For 
Qm(E) > 0, the L x L Green's matrix G^'m for 1 < n, m < N is defined by 

G%»{k,l) = (n,k\(H N -El)- l \m,l) , k,l = l,...,L. 
It follows from the Schur complement formula that 

GE,N f t t 7-t 1 rp f<E,N—l rp 

TV, TV — [ V N — & 1 — J TV^/V-l,TV-l J A 



-1 

■ N 
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Iteration of this formula gives a matricial continued fraction expansion: 

G%% = (v N -El-T N (---(V 2 -El-T 2 (V 1 -El + ()- 1 T 2 y 1 ---y 1 T N 

As Zf = E 1 — V\ — C, one sees that this is just the iteration of the Ricatti equation. Hence one 
deduces 

G E N ' N N = -(ZjJ)" 1 = -TT^)" 1 , 

where the inverse (Z^)" 1 exists because it is given by the Mobius transformation of Z^ G 
with ijE. SP(2L, C). Hence G^'^ G V^. Furthermore the geometric resolvent identity shows 

These identities allow us to note a few useful identities linking the entries of the transfer matrix 

N "' 1 ~ \ C D ) > 

to the Green's function: 

A-l _ r E,N n -\ _ r E,N-\rp r A-l _ n E,N j-ln _ r lE < N 

A - -Cr ljJV , U - -Lr lN _ l l N , U/i - —Lr NjN , AH — Lr^ . 

It is also possible to express B^ 1 , DB^ 1 and C _1 D in terms of Green's functions. 
3.7 Eigenvalue interlacing 

In this section, the above information on the spectrum of Jacobi matrices with matrix entries is 
complemented by a simple consequence of the min-max principle. In the case L — 1 of a Jacobi 
matrix, this is the theorem on alternation of zeros of the associated orthogonal polynomials. It 
also implies that the bottom (resp. top) of the spectrum of H N is less (resp. larger) than or 
equal the bottom (resp. top) of the spectrum of H N ~ l . 

Proposition 8 Let H N and H 1 ^" 1 be defined with Dirichlet boundary conditions on the right 
boundary. Then the eigenvalues of H N and H 1 ^ -1 satisfy the following interlacing property: 

E N < E N-1 < i = l,...,(iV-l)L. 

Proof. Let Hn = £ 2 (l, ■ ■ ■ , N) ® C l and the projection in Hn on the states on the right 
boundary, namely, in Dirac notation, on the span of (\N, /))«=i,...,l- Hence H-n-i — (1 — Hn)7~(-n 
and H-n-i C H-n with the natural embedding. Also f/ Ar_1 |-?/>) = H N \ip) for ip G H,n-i, i-e. 
UniP = (the natural embedding is suppressed in this notation). The min-max principle states: 

Ef = sup inf (V> | H N | V) , 

UcH N , dim(U)<j V , et/ ± , ||V!l=i 

where the supremum is over subspaces U of H,n, and the infimum over unit vectors in their 
orthogonal complement. For if^ -1 , the above facts imply 

£f _1 = sup inf | H N | V) , 

UcHn , dim([/)<7 , n N u=o ix^u 1 - , \\tp\\=l , U N i>=0 
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where the orthogonal complement is calculated in Hat. Hence the inequality Ej ~ > Ej follows 
because the condition UnU = is redundant and then one obtains a lower bound by dropping 
the constraint II v^ = 0. Next, for a subspace U C TCn with UnU — 0, let U = U © Hn'Hn- 
Then dim([7) = dim(C/) + L. Furthermore, the conditions if) G U 1 - and = are equivalent 
to ip E U L . Therefore, upon relaxing the constraints on the supremum: 

Ef- 1 > sup _ inf | H N | V>) = £f +L , 

which is the second inequality. □ 



4 Jacobi matrices with random matrix entries 

In this section we consider Jacobi matrices H N (u) with matrix entries uo = (V n ,T n ) n >i which 
are independent and identically distributed random variables drawn from a bounded ensemble 
(V a , T .) crg s of symmetric and positive real matrices. Expectation w.r.t. to their distribution will 
be denoted by E CT or simply by E. All formulas in Section I4TT1 also hold for more general covariant 
operator families and systems without time-reversal symmetry. Associated to each to are transfer 
matrices T^(u), Lagrangian planes $f (o>), their parametrizations Z^{oj) and U^(u), matrix 
entries A^(lj) and B^{uj) as in (f4"l~j) . total rotations, etc. In order not to overload notation, the 
index u is suppressed throughout. The basis change A4 will be taken independent of uj though. 



4.1 Integrated density of states and sum of Lyapunov exponents 

The integrated density of states of a random family of Jacobi matrices with matrix entries is 
defined bv [CSl IBLl ICL] 



M{E) 



lim — — 

N^oo NL 



E # { eigenvalues of H N < E } . 



According to ( I4"U1) . ( 142]) and the fact that C e = B e and D e = A e , one has 

1 N r E d 

M ^ = J im ¥7 E — 3m d e E log (det(^ + S« (U^T ) ) 



(46) 



For a fixed energy, this quantity can also be understood as a rotation number in the sense of 
Ruelle [Rue]. The second ergodic quantity considered here is the averaged sum of the positive 
Lyapunov exponents, denoted shortly by ^(E) here. For any complex energy E 6 C, it can be 
defined bv jUSl IBLl [ESI IUL1 ISE] 



l{E) 



1 



lim 



E log 



a 1 



n 



vn=l 



(47) 



where A L T is the L-fold exterior product (second quantization as for evolution group) of the 
symplectic matrix T, and the norm denotes the operator norm on the fermionic Fock space 
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A L C 2L . It is well known that j(E) is subharmonic in E \CS\ ICLj . Furthermore, the Thouless 
formula linking M and 7 holds [CSl IKS] . Actually this is the integrated version of the Kramers- 
Kronig relation stating that 'y(E) + iirN{E) for real E is the boundary value of a Herglotz 
function (which is given by the expectation value of the trace of the logarithm of the Weyl- 
Titchmarch matrix [KSJ, for which one has a matrix- valued Herglotz representation [GTj ). This 
is reflected by the following proposition showing that 7 and ttJV can be calculated as real and 
imaginary part of the Birkhoff sum associated to a single complex valued additive cocycle. 

Proposition 9 For E with Qm(E) > 0, 

1 - 

7 (£) = lim — ^ e E log ( det ( A% + B* (U^)*) ) . 



71=1 



Proof. Clearly one may replace in fj47|) by AiT^ because the boundary contributions 
drop out in the limit. Furthermore, instead of calculating the operator norm in (HT1) . one may 
insert a real Lagrangian plane $0 = (<j>i, ■ ■ ■ , 4>l) as initial conditions 



j(E) = lim — - E log 



a l \ T[mz? m- 1 ^ 



>i A . . . A L 



where now the norm is that of a vector in A L C 2i [BI4 A. III. 3. 4] (for covariant, but not necessarily 
random Jacobi matrices with matrix entries, this holds as long as E contains an average over <3> 
w.r.t. some continuous measure |JSS[ ISB] ). Recalling that the norm in A L C 2L is calculated with 
the determinant, it follows that 

j(E) = Hm ^Elog^det^*(nA4^ . 

Now ^n^=i^ ^rf-^ -1 ) ^0 = &n an d one ma y telescope (boundary terms vanish in the limit) 
and insert the Cayley transformation: 

7(£)= lim 4_ B f lor ( <•*««;>•<«?» 



00 2NL ^detftCe^Wi)), 

In each term, one can now apply Lemma E] for T = CMT^M^C* and $ = C$%_ v The 
hypothesis of the lemma are indeed satisfied for any E with ^sm(E) > because of the arguments 
in Section [2751 According to the definition of C/„ , it therefore follows 

1(E) = lim — - — E log ( dG f( {U / yU j + 1 \ I detUf + B* (17* J*) | 2 ) . 

The first contribution telescopes back again and the boundary term at N is bounded because 
1 < det(U*U + 1) < 2 L for every U G U 9^©^. Hence the first contribution vanishes in the 
limit. The second contribution is precisely the term appearing in the proposition. □ 
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4.2 Random perturbations 



This section gives the precise hypothesis of Theorem [2] and then provides the proof. Hence V n 
and T n are random and depend on a coupling parameter A > as described in the introduction 
and they give rise to transfer matrices T E {\) which depend analytically on A (lower regularity 
is actually sufficient). Throughout this section E G R. The first step of the analysis consists in 
the symplectic diagonalization of T E = T E {Qi) by an adequate symplectic basis change M.: 

M T E (X) M~ x = K exp(AP n + C(A 2 )) . (48) 

Here (and in matrix equations below) the expression 0(X 2 ) means that we have an operator norm 
estimate on the remainder. Furthermore V n is a (random) element of the Lie algebra sp(2L,R) 
calculated from (v n ,t n ) and 1Z is of the symplectic normal form of the free transfer matrix T E 
chosen as follows. The eigenvalues of T E form quadruples (A, 1/A, A, 1/A) which collapse to pairs 
(A, 1/A) if A G S 1 and A G R. If A G S , one speaks of an elliptic channel. Let there be L e of 
them, denote their eigenvalues by e 1711 , . . . , e tr,Le and set r\ = diag(r/i, . . . , t/lJ- As T e is supposed 
to be diagonalizable, the remaining Lh = L — L e channels are hyperbolic. The moduli of their 
eigenvalues are e Kl ,e~ Kl , with Ki > and for I = l,...,L h . Set k = diag(/«i, . . . , If a 

hyperbolic channel stems from a quadruple, it moreover contains a rotation by the phase of its 
eigenvalue A. This will be described by S G O(L^) which is a tridiaganol orthogonal matrix 
containing only either 1 or 2 x 2 rotation matrices on the diagonal and which satisfies [S, e K ] =0. 
The symplectic basis change M. is then chosen such that 



K = 



( Se K 








\ 





cos(t/) 





sin(r/) 








Se- K 







— sin (77) 





cos(?7) J 



Furthermore, let Ph and P e denote the projections (L x L matrices) onto the hyperbolic and 
elliptic channels. In particular, Ph + P e = 1 and diag(Ph, Ph) as well as diag(P e , P e ) commute 
with 1Z. The reader may consult |SB| where the basis change M. is constructed explicitly for the 
example of the Anderson model on a strip. Next let us state the precise hypothesis of Theorem [2j 

Hypothesis: The expansion factors K\ and rotation phases rji satisfy 

g h = min (l-e~ Kl ) > 0, g e = min |l - e l{m+m) I > 0. 

l<l<L h l<l,k<L e 1 1 



In order to develop the perturbation theory, some further notations are needed: 

r kAT*<\\ KA-^r* - ( A »( A ) B n( x ) \ _ ( A + \a n B + Xb n \ 2 
CMT n {X)M C - -^J - [ 1TTWn a + X^-J + ° {X 



Comparing with (1481) . one checks that 
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and 

«.= (s)' cp » c, (?) • 

Note that both A and P commute with both Ph and P e . 

According to Section 14.11 the averaged Lyapunov exponent and IDS at a real energy E for 
the A-dependent random operators H (A) are given by 

i * r E 

lx (E) + inJ\f x (E) = jim — - £ / de d e E log ( det ( A e n (\) + B°(X) (t^_ x (A))* ) ) . 

Let us expand the integrand 

log(det(Af(A) + Pf(A)(t/,f_ 1 (A)r)) = log(det(Af(A) + B E n {\)P h )) + J n (A) , (49) 
where 

J n (A) = Tr (log( 1 + (A E (\) + B E (X) P h )- l B E {\) (te,(A) - P h )* ) ) , (50) 

which is possible because A E (X) + B E (X) Ph = A + B + 0(\) is invertible for A sufficiently small. 
Note that only J„(A) depends on U E _ X , while the first contribution gives an contribution to IDS 
and Lyapunov exponent which can be readily calculated, similar as in (143]) . Hence we need to 
focus on the control of Birkhoff averages of J„ (A). For this purpose, one would first like to expand 
the logarithm in (j50|) . As B E {\) = BP h + 0(A), one therefore has to show that P h U E (X) - P h is 
small in norm. This means that the hyperbolic part of the dynamics (at A = 0) alines PhU E (X) 
deterministically with Ph up to small corrections due to the random perturbation. The following 
lemma is a strengthening of prior results [SBJ on this dynamical separation of hyperbolic and 
elliptic channels. 

Lemma 7 Let U E = 1 (choice of initial condition). Then there exist positive constants Ci,C2 
such that for A < -r^— and all n > 1, one has 

J 4ciC2 — 7 

\\P h u°(\)-P h \\ < 

9h 

Proof. For sake of notational simplicity, let P denote Ph within this proof. Let U be a symmetric 
unitary and set U' = CMT E {X)M.~ X C* ■ U . We will show that uniformly in n holds 

II PU' - P || < (1 - g h ) || PU - P || + ci A + c 2 \\PU-P || 2 . (51) 

An elementary dynamical argument then allows to conclude the proof (the constants C\ and c 2 
are then the same as in the statement of the lemma). In order to prove (15T1) . let us first note 
that A i — ► CJvlT E (\)M.~ 1 C* ■ U is an analytic path of unitaries. Hence the eigenvalues vary 
analytically in A |Kat j . Therefore U' = C1ZC* • U + R\ where Pi depends on U and T E (X), but 
one has the norm bound ||Pi|| < c\\ uniformly in n and U (because of the uniform bounds on 
the norms of V n and the error terms in (HHl) ). Therefore it is sufficient to show (ISTI) for A = 0. 
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By definition of the Mobius action, U' = (AU + B)(BU + Ay 1 so that 

PU' = P(A{PU-P) + (A + B))(A + B)- 1 (1 + B(PU -P) (A + B)~ l ) 1 . 

Now one can expand the last inverse in (PU — P) to first order, with an error term bounded by 
1 1 PU — P\\ 2 . Then multiplying out all the remaining factors shows that 

PU' -P = P(A-B) (PU -P)(A + B)- 1 + R 2 , (52) 

with an error term that satisfies ||i? 2 || < c 2 \\PU-P\\ 2 . Now \\(A + B)- 1 ]] = 1 and \\P(A-B)\\ = 
maxi<K£ h e~ K ' = 1 — which implies (}5Tl) for A = 0. □ 

Now it is possible to expand the logarithm in (!50|) because B^(X)P e = 0(A) so that Lemma[7] 
implies B^(X)(U^_ 1 (X) — Ph)* = 0(X) (here still all error terms are norm bounded). Hence 

J n (X) = Tr((Af(A) + B^(X)P h )-'B^(X)(Ul 1 (X)-P h y) + o(—) , 

where the L comes from carrying out the trace after having applied the norm bound. Expanding 
A^(X) and B^(X), using the commutativity of A and B with P h and P e and invoking Lemma [7] 
in order to show PhU^(X)P e = 0(X/gh) now implies 

LA 2 ' 



J n (X) = Ti((A + B)- l {(P h Uti(^TPh-P h ) + Xb n P e Uti(\yPe)) + O(-) 
v J v 9h ' 

Setting 



N-l N-l 

h(N) =E-—J2 {PhU*(\)P h - Ph) , h(N) = E — £ PeU*(X)P e 



n=0 n=0 



one has 

JV 



E^ J "( A ) = Tr ( (A + B)~ 1 Ih(N)* ) + ATr((A + J B)- 1 E (T (6 (T )/ e (iV)*) + 



In order to calculate and bound the two traces, we will use the Hilbert-Schmidt spaces Tih 
and 7i e of complex matrices respectively of size x and L e x L e , furnished with the scalar 
product (C\D) 2 = Tr(C*D). The corresponding norms are ||C|| 2 = Tt(C*C)2. They satisfy 
norm the inequality ||C||2 < yL ||C|| w.r.t. to the operator norm (where one may, of course, 
also use respectively Lh and L e instead of L). For a L x L matrix C, we will identify PhCPh and 
P e CP e with vectors in respectively TCh and H e . Let us first focus on the second trace in the last 
expression. The Cauchy-Schwarz inequality implies 

TrdA + B^E^IeiN)*) < VZWP^A + By'E^PeW \\I e (N)\\ 2 . 

As the operator norm appearing on the r.h.s. is bounded, the following lemma shows that this 
trace is of order 0(X/y/L) and hence does not contribute to leading order. 
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Lemma 8 There exist positive constants c\ and C2 such that 

IWiV)ll2£ ^( ciA + 4 

Proof. This is a matrix version of the oscillatory sum argument in |PFl ISBj. First note that for 
each summand in I e (N), one has P e U^(X)P e = e iri P e ^_ 1 (A)P e e 1 ' 7 + (9(A). Thus 

I e (N) = e^I e (N)e^ + R x + R 2 , 

with an average error term Ri satisfying ||i?i|| < C\\/L and boundary terms R2 satisfying 
\\R2W < c 2 /(NL). Now let us define the super-operator D v : H e -> H e by D V (C) = e lr? Ce lr? . This 
operator is diagonal and the hypothesis g e > implies that (1 — D v )~ l exists and its norm is 
bounded by l/g e . As I e (N) = (1-D V )- 1 (R 1 +R 2 ), it foUows that ||ie(A0|| 2 < (H-Ri II2 + \\Rah)/g e 
which leads to the desired bound. □ 

A similar argument allows to calculate the remaining trace. 
Lemma 9 One has 

TidA + B^hiNy) = 2X^eTi((e 2K -l)- 1 P h E a (a a + b a )P h ) + , -L^ . 

Proof. One first has to refine ( 1521) and include the O(X) contribution. Invoking Lemma [7] at 
several reprises, some lengthy but straightforward algebra shows 

P h U*(X)P h -P h = Se~ K {P h Utx{\)Ph-Ph) e' K S f + 2X^eP h E a (a a + b a ) P h e' K S f + LR lt 

with H-Rill < C\X 2 j{Lg\) and the formula is understood as identity for operators on Ti^. Now 
define the super-operator D K : Hh — > 'Hh by D K {C) = Se~ K Ce~ K S t . One directly checks ||(1 — 
Ac) -1 !! < As above, 

I h (N) = (1 - DJ- 1 ^ Vie P h E a (a a + b a ) P h e- K S f + R x + i? 2 ) , 

where H-R2II < oi/{NL). Replacing this into the trace and bounding the error terms by the 
Cauchy-Schwarz inequality just as before Lemma M allows to bound the error terms. The leading 
order contribution can be calculated using the identity 

S l e- K (1 -D*)- 1 (Se~ K ) = (e 2K - l)" 1 . 

This completes the proof. □ 
Replacing Lemma [9j it then follows that 

lx {E) + nrNx{E) = E a \ f E ded e log(det(^(A) + B e a {X) P h )) 

L J-00 

(53) 

+ ^ sfte Tr ( (e 2K - l)~ 1 P fc E a (a CT + K) P h ) + O (- , ^ V 

\9e g h j 
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As this holds also for the translation invariant operator with one fixed a at every site n (which 
has non-random independent entries), one deduces 




where ^\^(E) and Mx,a{E) are averaged Lyapunov exponent and IDS of the translation invariant 
operator (as already described in the introduction). This therefore proves Theorem [2J Expanding 
the logarithm in (153"]) . some straightforward algebra leads to more explicit perturbative formulas 
for 7a(-E") and N\(E) in terms of E a (V a ). Let us note that the contribution of Lemma [9] is real 
and hence only contributes to / ~fx{E). The lowest order contribution to Afx(E) is given by the 
first term in fj49]) and one can check that it is only the contribution of E(P .) which changes the 
rotation phases of the elliptic channels of TZ. 
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